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Abstract 

In previous work the authors analysed the global properties of an ap- 
proximate model of radiation damping for charged particles. This work 
is put into context and related to the original motivation of understand- 
ing approximations used in the study of gravitational radiation damping. 
It is examined to what extent the results obtained previously depend on 
the particular model chosen. Comparisons are made with other models 
for gravitational and electromagnetic fields. The relation of the kinetic 
model for which theorems were proved to certain many-particle models 
with radiation damping is exhibited. 



1 Introduction 

The study of gravitational radiation in general relativity relies in an essential 
way on sophisticated approximation methods. It is important to understand the 
relation between the approximate models and the exact theory. Many of the 
issues arising in this context have recently been surveyed by Blanchet The 
relation between the exact and approximate equations is relatively straightfor- 
ward. What is much more difficult is to establish a relation between solutions 
of the two sets of equations. 

A standard approach to the theory of gravitational radiation (see e.g. 
' d] ) is to combine post-Newtonian and post-Minkowskian approximations via 
matching. There are some mathematical results available on these individual 
types of approximation. It has been shown that under suitable conditions the 
post-Minkowskian approximations are asymptotic |^ but divergent ]l9| . The 
post-Newtonian approximations are more difficult to handle. A general mathe- 
matical framework for formulating the relevant questions was set up in pO|. It 



was shown how the weh-known divergences of higher order post-Newtonian ap- 
proximations could be understood within that approach. This framework also 
provided a basis for showing that the post-Newtonian expansion is asymptotic 
at the lowest order (Newtonian) level It is reasonable to expect that anal- 
ogous results could be proved for the IPN and 2PN levels but this has not been 
done yet. 

None of the results mentioned above prove anything about matching the 
two approximations and until this can be done the possibilities of understand- 
ing anything about radiation (even at the quadrupole level) in a rigorous way 
are very limited. In the following these limits will not be removed. A more 
modest goal is pursued, which is to understand some of the effective models 
for radiation themselves while postponing the question of their relation to the 
exact equations. Radiation comes in at the 2.5PN level, but there are hybrid 
models where radiation at that level is combined with a description of matter 
and non-radiative gravitational fields at the Newtonian or IPN level. These 
models will be referred to as (0-f 2.5)PN and (1-1-2. 5)PN models. 

For purposes of comparison it is valuable to consider the case of charged 
matter emitting electromagnetic radiation as an analogue of matter emitting 
gravitational radiation. In this case the post-Minkowskian expansion becomes 
trivial due to the linearity of the Maxwell field while there is a non-trivial ana- 
logue of the post-Newtonian expansion. In the electromagnetic case radiation 
first appears at the 1.5PN level and it is possible to define (0-1-1. 5)PN and 
(1-1-1. 5)PN models. These will be referred to as the Coulomb and Darwin cases 
respectively. 

In ||l5| the authors proved theorems on the global existence and asymptotic 
behaviour of solutions of a certain model of the (0+1.5)PN type in the electro- 
magnetic case. The purpose of the present paper is to put this mathematical 
result into a wider context and to discuss its ramifications. For technical details 
of the proof the reader is referred to [|l^ . 

A general mathematical setting for expansions of PN type is as follows. Mat- 
ter is described by a function /(A) depending on a parameter A belonging to 
an interval (0, Aq]. The field (gravitational or electromagnetic) is described by 
a function 5(A). The idea is to consider a one-parameter family of solutions 
which become non-relativistic in the limit A — > in some appropriate sense. 
The parameter A corresponds to where c is the speed of light. More pre- 
cisely, the given functions / and g describe a one-parameter family of physical 
systems represented in suitable parameter-dependent units. Then c is to be in- 
terpreted as the numerical value of the speed of light in the given units. A more 
conventional physical description of the post-Newtonian expansion would say, 
for instance, that velocities of the matter are small with respect to the speed 
of light. The connection between the two descriptions is that to say that the 
matter velocity is of order unity and the speed of light large in one system of 
units is equivalent to saying that the matter velocity is small and the speed of 
light has a fixed value in a different (fixed) system of units. 

In our discussion of expansions for charged matter we will start with the 
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following formal expressions: 

E = Ea + X^'^E^ + \E2 + X^'^E:^ + . . . , (1) 
B = Bo + X^/^Bi + \B2 + A^/^^a + . . . , (2) 
P = Pa + A'/'pi + Ap2 + A^/^s + • • • , (3) 

j = JO + Al/271 + Aj2 + A3/2j3 + . . . . (4) 

The use of expansions in half-integral powers of a parameter A ensures that the 
notation fits with the usual terminology of post-Newtonian expansions. If we 
write A = then the Maxwell equations are: 

div£; = 47rp, curli; = -X^^^dtB, (5) 
divB = 0, curls = X^^^dtE + inX^/'^j. (6) 

Substituting the expansions into the equations and comparing coefficients gives 
a sequence of equations for these coefficients. Note that divBg = and curli^o = 
0. It follows that assuming fields defined on R'^ and vanishing at infinity Bq = 0. 
This condition will be assumed from now on. The first equations for coefficients 
of E of the sequence (ordering by powers of A) are divi?o = 47rpo and curli?o = 
0. The second equation implies that Eq is the gradient of a function U and 
then the first equation implies that AU = AirpQ. Thus the Poisson equation of 
electrostatics is recovered, with U being the electrostatic potential. 

Next a remark will be made concerning the coefficients of the odd half-integer 
powers in the expansion of E, p and j and the integer powers in the expansion of 
B. These can be consistently set to zero. In fact if the coefficients belonging to 
this class are set to zero in the matter quantities then (assuming fields defined on 
R'^ and vanishing at infinity) the coefficients of the fields belonging to this class 
automatically vanish. In this section it will be assumed that all these coefficients 
are identically zero. They will, however, be important in the later section |^. In 
a pure post-Newtonian expansion like the above they play no essential role and 
if they are maintained they just lead to 'shadow equations' as discussed in |^ . 

Consider next Bi. It satisfies the equations divi?i = and curlBi = dtEo + 
'in jo. Combining these equations and using the fact that curli?o — divBi = 
gives ABi = — 47rcurljo- Thus Bi solves a Poisson equation with a source deter- 
mined by the lowest order coefficient in the expansion of the current. It is easily 
shown, using the asymptotic conditions assumed above, that the solution of this 
Poisson equation satisfies the original first order equations. Next combining the 
equations divi?2 — 47rp2 and curl_B2 = —dtBi and substituting for curlBi leads 
to the equation 

AE2 = 92^0 + 4^(Vp2 + dtjo). (7) 

Assuming that the matter variables have compact support for each fixed value 
of t the second term on the right hand side of this equation also has compact 
support and leads to no difficulties in solving this Poisson equation. It remains 
to discuss the first term, which is not compactly supported. Assume that the 
matter sources satisfy the continuity equation 9tp + divj = 0. Then the total 
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charge, / p{x)dx, is independent of time and the first term on the right hand 
side of the Poisson equation for E2 is 0{r~^). It follows (for example from 
results in the appendix of |2^]) that there is a solution for E2 which vanishes 
at infinity and is unique in a suitable function space. The case of B3 is more 
problematic. According to the assumptions made up to now it should satisfy 
the Poisson equation: 

AB3 = d^Bi - 47rcurlj2. (8) 

Once again the second term is unproblematic, but the first needs to be looked 
at more closely since it might have a contribution which only falls off like 
as r 00. This contribution is made by the second time derivative of the 
quantity J jo{x)dx. Using the continuity equation and partial integration the 
latter quantity can be rewritten as —{d/dt) J xpo{x)dx. This is (up to sign) the 
first time derivative of the dipole moment D{t) of the charge distribution at the 
Coulomb level. Thus the coefficient which is problematic is d^D/dt^ and this is 
a quantity which comes up in the description of radiation damping. What has 
happened here, roughly speaking, is that with B3 which is the coefficient in a 
term of order \^^^ we have reached the level at which radiation damping plays 
a role and the naive post-Newtonian expansion breaks down. A correct analysis 
using matching would show that the asymptotic conditions to be imposed on 
i?3 are different from the vanishing of that quantity at infinity. With a right 
hand side of the given form the Poisson integral will diverge and it is not to be 
expected that the Poisson equation for ^3 will have a solution which vanishes 
at infinity. In any case it does not have a solution in the kind of weighted spaces 
introduced in ||2^. Thus the expansion here breaks down at the 1.5PN level in 
much the same way as the post-Newtonian expansion for gravity breaks down 
at the 3PN level (as shown in pO|). 

Adopting terminology adapted to that used in we call a solution of the 
equations for £^0 the Coulomb approximation, a solution of the equations for Eq 
and Bi the quasi-electrostatic approximation, and a solution of the equations 
for Eq, Bi and E2 the Darwin approximation. 

In order to throw some more light on the particular form of the expansions 
chosen it is useful to reexpress the Maxwell equations in four-dimensional form. 
The Maxwell equations in Minkowski space are 9a F"^ = Airjp and daFp^ + 
+ dpF-^a — 0. Here F"'^ is the field tensor and j" is the four-current. 
The expansions above can be written in terms of one-parameter families of 
these quantities. Since we are interested in the Newtonian limit it is useful to 
introduce a family of metrics /a/3(A) as in by /oo = -1, /oa = 0, Jab = ^5ab- 
All these metrics are flat and in fact represent the Minkowski metric in different 
units. As basic quantities we take one-parameter families F"^(A) and j"(A) 
defined for A > and assume that they have limits for A ^ 0. The tensor 
/a/3 (A) also has a limit for A — > but the limit is not a Lorentz metric. In 
general indices will be raised and lowered using the metric fap and its inverse. 

Define E"" = F"o, B^ = \-^/^F^3, B^ = X-^/^F\, B^ = X~^/^F\. The 
last three relations can also be expressed in terms of the spatial volume form 
associated to /, call it e/. Using this we can write = (e/)°f,^F^c. In general 
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situations will be considered where F°o has a non-zero limit and then the same 
is true of E"". On the other hand it will be assumed that F°'i, is 0(A) as A — > 
and then 5° is automatically 0{\^^^). Define p — and take the current in the 
three-dimensional form of the Maxwell equations to have components With 
these definitions the three- and four-dimensional forms of the Maxwell equations 
are equivalent and the expansions of E, B and j above correspond to expanding 
F°'fS in integral powers of A, starting with the power zero, and assuming that 
the first coefficient in the expansion is zero when both a and /3 are spatial. The 
correctness of this statement depends on the fact that the coefficients of the 
odd half-integer powers in the expansion of E, p and j and the integer powers 
in the expansion of B have been set to zero. Radiation damping comes in at 
order A'^/^ and this is not allowed for by the expansion in integral powers of A. 
However as long as we are doing an expansion of post-Newtonian type it is the 
case that if j" contains only integral powers of A the same is true of i^"^. 

In the practical applications of post-Newtonian approximations it is usual to 
expand the fields but not the matter quantities. This point has been discussed 
briefly in section 6 of (20j . What this means in effect is that we are not talking 
about functions which solve a certain system of equations exactly but rather 
about parameter-dependent families of functions which satisfy certain equations 
up to an error of the order of a power of the parameter. Related to this is the 
fact that the equations themselves are often only fixed up to an error of a certain 
order. In reality we are dealing with equivalence classes of equations. This point 
will be illustrated by comparing the Darwin model introduced by Degond and 
Raviart in pO| with the Darwin approximation introduced above. Degond and 
Raviart split the electric field into transverse and longitudinal parts and set 
the transverse part of E to zero in the Maxwell equation containing curlB to 
obtain the Darwin model. We claim that at the level of what we call the Darwin 
approximation the expansions of the Darwin model of |10 and the full Maxwell 
equations are identical. Thus at that level both sets of equations are in the 
same equivalence class. In terms of our expansion the transverse part of E is 
that arising from the equations for cwc\E„ while the equation containing curlB 
gives rise to the equations for curli3„. The only coefficient Bn which occurs at 
the level of the Darwin approximation is Bi , in the equation curli?i the only 
contribution from the electric field comes from Eq, and Eq has no transverse 
part. Thus on this level the passage from the full Maxwell equations to the 
Darwin model has no effect. 

Up to this point we have concentrated on the fields and said little about the 
matter. The equations of motion for charged particles contain the Lorentz force 
Fl = E+\^^^v X B. In terms of the expansion coefficients this gives i?o + A(i?2 + 
V X Bi) as the contribution up to order A. The most common choices of matter 
models when discussing approximation schemes are point particles or a perfect 
fluid. In general relativity the concept of a point particle is very problematic 
and can itself only be defined in an approximate sense. A perfect fluid is better, 
but when it comes to obtaining mathematically rigorous results it also has its 
problems. On the one hand there is the problem of the formation of shocks, 
which is always a danger when considering solutions of the Euler equations on 
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a long time scale. When dealing with radiation we do wish to consider the 
evolution of the system on a long time scale. Once shocks are formed we enter a 
regime where little control is possible with known mathematical techniques. The 
other problem concerns the question of fluid bodies. In studying radiation we 
want to consider isolated systems. The easiest thing would be to take the fluid to 
have compact support. Unfortunately the Cauchy problem for solutions of the 
Euler equations with compactly supported initial density is poorly understood, 
even locally in time. There has been some progress, since Wu proved a local 
existence theorem for the incompressible Euler equations with free boundary 
in the irrotational case. Christodoulou and Lindblad proved estimates for 
the case with rotation which will hopefully lead to an analogous result in that 
situation. However much remains to be understood in that area. It might be 
thought that replacing compactly supported initial data with data which fall off 
at infinity might simplify the Cauchy problem for a fluid body but nobody has 
managed to take advantage of this yet. 

In Section 2 models of radiation which combine expansions at PN levels 
which belong to different powers of the expansion parameter ('hybrid models') 
are discussed. After reviewing the results of on one model of this type for 
electromagnetic radiation we discuss various issues relating to the possibility of 
generalizing those results to related models for electromagnetism and gravita- 
tion. In the model of the matter is described by kinetic theory. Section 3 
of this paper relates the kinetic model of to the underlying dynamics of a 
many-particle system with radiation reaction, thus providing a better justifica- 
tion for the kinetic model which was originally introduced in an ad hoc manner. 
The paper concludes with a brief outlook. 



2 Hybrid models 

By a hybrid model we mean a model where an expansion of post-Newtonian type 
up to a certain order is combined with effects of radiation which are formally 
of higher order in the expansion parameter A. Some variants of this procedure 
were already mentioned in the introduction. The model studied in ||l5[| was for 
the case of the electromagnetic field with a kinetic description of the matter at 
the Newtonian level and dipole radiation at the 1.5PN level. It was inspired by a 
model of Blanchet, Damour and Schafer [|| for the case of the gravitational field 
with a hydrodynamic description of matter at either the Newtonian or IPN level 
and quadrupole radiation at the 2.5PN level. The replacement of the fluid by 
kinetic theory was motivated by the difficulties with a mathematical treatment 
of fluids mentioned in the introduction. The replacement of the gravitational 
by the electromagnetic field was motivated by simplicity. The extent to which 
the particular choice of model in [ p^ was essential in obtaining results will be 
discussed below. 

In [Q an initial model with higher order (in fact fifth order) time derivatives 
which could not be expected to have a good initial value formulation was modi- 
fied by two methods for reducing the order of the time derivatives. The resulting 
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model is not equivalent to the original one but is in the same equivalence class 
in a sense already indicated. In other words the equations of the two models 
agree up to terms which are regarded as being of higher order at the given level 
of approximation. One method of reduction is to do an appropriate change of 
matter variables. The other is to substitute the equations of motion into the 
undesirable time derivatives and to discard terms which arise from this which 
are formally of higher order. In Q the first method is used to remove two time 
derivatives and the second to remove the remaining three. Thus we could call 
this a (0+2.5)PN or (H-2.5)PN model with 3+2 reduction of time derivatives. 
The model of |l5) is a (0+1.5)PN model with 2+1 reduction of time derivatives. 

The results of |]l5| will now be stated. For reasons explained in that paper 
it is necessary to have species of particles with different charge to mass ratios in 
order to get an interesting result. The particular choice made was to have two 
species of particles with unit mass and charges of unit magnitude and opposite 
sign. The phase space densities of these two species of particles are described 
by functions /+ and /~. These functions satisfy the Vlasov equation with a 
force term which is the sum of a Coulomb term resulting from a potential U 
and a radiation reaction term. The primitive form of the latter term is (up to 
sign) e(PD/dt^ where 

D{t) ^ J x{p+{t,x) - p^{t,x))dx (9) 

is the dipole moment of the charge distribution and e is a small constant. The 
charge distribution itself is defined to be 



p{t,x) ^ p'^{t,x) - p {t,x) = J r{t,x,v) - f {t,x,v)dv. (10) 

The primitive form is reduced by procedures of the type already mentioned in 
such a way that the equations for and are of the form 

dtf+ + {v + ei^I'l (t)) • V,/+ + WU ■ V,,f+ = 0, (11) 

dtr + {v- (t)) • V,/- - VU ■ V,/- = 0, (12) 

where 

D^^^{t)^ I \/U{t,x){p+{t,x) + p-{t,x))dx. (13) 



The full system of equations for U and / is called the VPD system (Vlasov- 
Poisson with damping). Define the total energy £ of the system by 



£(i) =fki„(i)+^pot(t), (14) 



with 



fkin(<) = y \'v\'^{f^it,x,v) + f {t,x,v))dxdv, and 

fpot(t) = \'^U{t,x)\Ux, (15) 
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denoting kinetic and potential energy, respectively. With these preliminaries in 
hand we can state the main theorem of | p5| . 

Theorem 1 // are smooth initial data with compact support for the VPD 
system, then there is a unique smooth solution of VPD for t > with data 
f^(t = 0) = f^. The energy evolves according to 

£it) = -e\Di'\t)\\ (16) 
Moreover, the following estimates hold for t G [0, oo[. 

(a) < C{1 + /or p G [1, §]; 

(b) ||VC/(i)||p<C(l + t)-^ /orpG [2,f]; 

(c) <C{l + t)-^. 
Here \\ \\p denotes the norm. 

Thus we obtain a global in time existence theorem and statements about how 
various quantities decay as t oo. The theorem was proved by adapting 
methods previously applied to the usual Vlasov-Poisson system by Lions and 
Perthame ||l8| (global existence) and lUner and Rein (asymptotic decay). 

Next we consider whether similar results can be obtained for other hybrid 
models. First we give the equations for some of these models. Starting with the 
primitive form of the equations for charged particles described by the Vlasov 
equation at the (0+1.5)PN level, a 1+2 reduction can be carried out. This means 
introducing new variables in the characteristic equations hy — X^ ^ ^D{t) 
and — =F ^D{t) with the result that the new characteristic equations are 

X = F± and V = -VC/(t,X± ± el)(t)). The Vlasov equation with these 
characteristics has the form 

dtf^ + V ■ V,/± - VC/(t, X ± eDit)) ■ V„/± = 0. (17) 

Then the other method must be applied to replace D{t) by D^^\t), where 

/ / v{f+ - f-)dxdv, 



cf. section 1 of ||T5[. This defines the model. There is no obvious obstacle to 
proving local existence for this case. The other hybrid model obtained from the 
same starting point is that resulting from a 3+0 reduction. In that case the 
Vlasov equation is of the form 

dtf^ + V ■ V,/± ± (Vt/(i, x) + e7?[3] (i)) . v,/± = (18) 

where 

D^^\t)^ I I VU{f+ + f-)dxdv+ I I VVU ■v{f+ + f-)dxdv. (19) 
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In this case proving local existence may be more difficult, due to the second 
derivatives of U occurring in the Vlasov equation. There is a loss of one deriva- 
tive compared to the 2+1 and 1+2 reductions. 

The consequences of the lost derivative will now be examined in some more 
detail. In order to get a local existence theorem the orders of differentiability of 
the different unknowns must fit together properly. Suppose that t/ is fc times 
differentiable. Then /, and hence p will be fc — 2 times differentiable. Thus, in 
order that a match be obtained, two derivatives must be gained when solving 
the Poisson equation. It is well known (see e.g. |l^) that this is not true in 
the context of pointwise differentiability properties. It is true when regularity is 
measured in Holder or Sobolev spaces |l^ . Thus in order to get a local existence 
theorem it is cither necessary to work with spaces of that kind or to see that at 
some point a derivative can be gained in comparison to the naive counting just 
given. This appears as a borderline case for local existence. 

In the case of the 2+1 reduction studied in a simple and convenient for- 
mula for the energy loss of the system due to emission of radiation was obtained. 
For the 3+0 and 1+2 reductions things seem to be more complicated. In the 
model considered in |]l5| the rate of change of the total energy is —e\D^'^\t)\^ , 
which is manifestly negative. The corresponding expression for the 3+0 re- 
duction is tD^^\t) ■ D^^\t), which is neither manifestly negative nor obviously 
related to the other expression. On a formal level one link can be seen. Suppose 
that we replace the quantities Z3t"l(i) by the true derivatives D'^"'\t). Then we 
have the identity that 

D^^\t) ■ D^^\t) = -\D^^\t)\'^ + dt{D'^^\t) ■ D^^\t)). (20) 

If we suppose that the function D{t) has a suitable oscillatory behaviour then 
the last term will average to zero so that in an average sense the left hand side 
is equal to the first term on the right hand side. To turn this formal argument 
into a rigorous one would require at the very least a lot more information than 
was obtained in jl^. The relation to the energy loss in the 1+2 reduction is 
even less clear. 

Consider next a model with gravitation and kinetic theory at the (0+2.5)PN 
level. In comparison to the electromagnetic case the main difference is that the 
third derivative of the dipole moment is replaced by the fifth derivative of the 
quadrupole moment. The primitive form of the radiation reaction term is (up 
to sign) e{(f'Qab/ dt^)x°' , where the quadrupole moment is given by 

Qab{t) = j [xaXb - \x\^5abj p{t: x) dx. (21) 

The fifth order derivative is reduced to the third order derivative by the trans- 
formation (cf. equation (10) in p5||) 

X ^ X - eQ(^\t) ■ X, (22) 
V^V- eQ(4) {t)-X + eQ(3) . y_ (23) 
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The transformed Vlasov equation is 

dtf + {v-2eQ^^\t)-v)-V.J 
-[V[/(t,a; + eQ(3)(t)) - eQ^^^t) ■ yU{t,x)] ■ V„/ = 0. (24) 

In doing this transformation terms which arc formaUy of order have been 
discarded. The next step in the 3+2 reduction of the gravitational case is to 
replace the time derivative Q^^^ by the reduced quantity Q^^l. This gives rise to 
second derivatives of the potential as in the 3+0 reduction of the electromagnetic 
case. To write the expression for the reduced time derivative Q^^^ in a compact 
form it is useful to introduce the following operation on tensors 

STFiTab) = (l/2)(Tafc + Tba) - il/3)6^''T,d5ab- (25) 
It takes the symmetric trace-free part of a given tensor. Then 

Q[^1(<) = STF^2y j[{v -VVU ^VU)®x + ZVU ®v\f{t,x,v)dxdv^. (26) 

This may be compared with equation (6.6) of Q| where a corresponding equation 
is given for the hydrodynamic case. Two of the terms above can easily be seen 
to correspond to terms in the equation of Q . The sum of the remaining terms in 
the latter equation correspond to the term above involving U via the equations of 
motion. There does not seem to be a useful way to eliminate the time derivative 
of the potential in the case where the matter is described by kinetic theory. 

As in the 3+0 reduction in the electromagnetic case the second derivatives 
may make proving local existence difficult. Note that a 4+1 reduction would 
lead to the appearance of third derivatives of the potential, the loss of two 
derivatives compared the model treated in jist and, presumably, to ill-posedness 
of the evolution equations. 

3 Comparison with A^-particle systems 

The VPD model considered above has been introduced as an effective model for 
electromagnetic radiation. The purpose of this section is to see how this model 
relates to the effective equations for TV-particle systems of Abraham-Lorentz 
type ||l^, |l^, |2^ in the limit oo. As before we assume the fields E and B 

do satisfy the Maxwell equations (|^), and they are expanded as in (^, (^ 
in A^/^ = . Contrary to section^ however, we do not suppose a priori that 
the coefficients of the odd half-integer powers in the expansion of p and j 
and the integer powers in the expansion of B are set to zero. 

3.1 The Darwin approximation 

In section |l| we have already seen that the first few equations the coefficients 
have to satisfy are i?o = and 

Ai?o = 47rV/9o, ABi = -47rcuryo, 

l\Ei = 47rV/9i, A^Ea = AiiV p2 + dt{dtEo + 47rjo). 
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We intend to derive an explicit approximation of the Lorentz force Fl = E + 
A^/^w X B up to order 0{X) — 0{c~^), given through 

Fl=Eo + X^'^Ei + X{E2 +vxBi) + 0{X^^^). (27) 

Suppressing the time argument, clearly 

Eo{x) = 47tA-\S/po){x) = / po{y) dy, (28) 

J \x-y\-^ 

with the analogous formula for Ei by replacing po by pi. Moreover, 

Bi{x) = ~47rA-\cxirlj„){x) = - / x jo{y) dy 

J \x~y\-^ 



x-y 
\x - yY 



xjiy)dy + 0{X'/^). (29) 



To determine E2, we note that dtpo + divjg = in view of divi?o — 47rpo and 
curli?i — dtEo + 47rjo, whence 

AidtEo) = 47rV(atpo) = -47rV(divjo). 
Accordingly, E2 satisfies 

E2 = 47rA-i(Vp2) + 9t A~'(at£;o + 47rjo) 

= 47rA-i(Vp2) + 47rat ( - A-2(vdivjo) + A-^jo 



Here we follow the device of |Q, p. 296] how to solve iterated Poisson equations, 
and observing that X{x) = J \x — y\ a{y) d^y is the solution to A'^X = — Sttct, 
after some calculation the more explicit form 

E2{x) = AnA-\Vp2){x)-^dt(^l j^^My)dy 

^—^ix-y)-jo{y)ix-y)dyj (30) 



\x - y\ 

is found. Conversely, it may be verified directly that this function E2 obeys 
divE2 = 47r/92 and curl£;2 = -dtBi. Taking into accoun t (p8| ), (||), and (|^, 
we have obtained an explicit representation for F^ from (^7|) up to order 0{X). 

To make the connection with the N 00 limit of a system of N particles 
we need to rewrite E2 further and make use of the fact that this description is 
not unique since terms of order ©(A"^/^) can be added or subtracted. Adopting 
the terminology from section |l| this means that we look for a more convenient 
representative in the same equivalence class of equations. From now on we 
assume that f = fo + X^^^fi + A/2 + . . . satisfies the relativistic Vlasov equation 

dtf + v-V,f + {E + X^/^vxB)-Vpf^O (31) 
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coupled to the Maxwell equations (||) , (|^) . Here p G R'^ denotes the momentum, 
thus V = {1 + \p^)~^^^p for the velocity if we consider particles of unit mass. 
Therefore v — p + 0{X) implies that 

= at/o + P • V,/o + ■ Vp/o + 0(Ai/2) (32) 

i.e., to lowest order / = /o is a solution of the usual Vlasov-Poisson system 
with electric field Eq obeying /S.Eq — AirV pq. Hence we can invoke (|3^) and 
jo = / w/o dp to write 



dth = - I P{P- V,/o + Eo ■ Vp/o) dp + 0{X^/^) 

p{p-Vccfo)dp+ f Eofodp + 0{X'^^). (33) 



Next we explicitly perform the c^j-differentiation in ( pOj ) and substitute the ap- 
proximate expression for dtja from (|33|). We obtain after some simplification, 
including Ea ^ AnA-^Vp) + 0{\^/^) and fo = .f + 0{X^/^), that 

E2{x) = AttA-\Vp2){x)-2ti [ [ ^^^-\Vp){y)J{y,p)dydp 



\x - y\ 

'^^ I i T^^{x-y)-A-\Vp){y)fiy,p)dydp 



\x - y\^ 

I f f f p" 3ip-{x-y) f 



2jJ\\x-y\^ |.-# ji--y)fiy^p)dydp 

+ C>(Ai/2). (34) 



In (prh the field E2 carries the factor A. Whence we can as well use E2 in its 



form ( p4[ ) rather tha n (|30| ) to obtain a valid approximation of Fl up to order 
0{X). In particular, (|34|) .s more convenient to compare this approximation to 
what comes out from the associated particle model. In it was show that 
the dynamics of N particles coupled to their self-generated Maxwell field can 
be described over long times by means of the effective equations 

^Tki„(w„) = GS°(g,«,«) + 0(A3/2), l<a<N, (35) 

where T^in is the kinetic energy, and the force G^^ is given by 

+ Ae„ 2^ e;3 - — - V(3 - 

Hvi3 ■ ^c^pf ^ K ■ ^ , (^g ■ gg/j) ^. ^ 
13 Up 15 - \p 13 ^^Z' + „I3 I • (36) 
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Here particle a has position qa{t) G and velocity Va{t) — qa{t)- Moreover, 
= Qa—qp- Note that in Lemma 3.5], eq. (^5|) is formulated on a different 
scale by means of a dimensionless parameter e ^ 1. In this units, the error was 
©(e''/^), and to pass to ( p5| ) it is first necessary to multiply by to undo 
the transformation, and then e = A, whence the error term is of order ©(A"^/^). 
In addition, in the charge distribution p and the current j in the Maxwell 
equations are considered without the factor in. Hence to adjust to (j^) 
the expression for G^^ had to be multiplied by iir. We also remark that (p5 
holds in the limit of slow particles which are far apart, cf. [p^ . The equation 
■^Tkin{va) = {q,v,v) is governed by a suitable Lagrangian, consisting of a 
Coulomb and a Darwin contribution, cf. Ch. 12.6] and |l^, Eq. (1.1)]. Since 
the Abraham model considered in |l^ is only semi-relativistic, the particular 
form of Tkin cannot be expected to agree with the corresponding expression 
resulting from the Vlasov-Maxwell system through expansion in powers of A^^^; 
see 1^ for a recent improved model. However, in the formal limit N ^ oo the 
force coincides with Fl = Eq + \^''^Ei + \{E2 + v x Bi), using i?2 from 
(^4[). To see this, we fix one particle a, i.e., its position x = qa € R'^ and its 
velocity v ^ Va G R"^. Setting — 1 for all /3 corresponding to (^, then 

results in the Coulomb part; note that the corresponding expression in is 
given by Eq + X^'^Ei + A47rA-i(Vp2) = 47rA-i(Vp) + 0{\^/^). In addition, 
we obtain 



/3=1 

I3=£a 



cf. (g9|). To deal with the remaining part of ( pq ) we have to reexpress the vp 
through lower order terms. Since A(— 2 | / g | ~ ^"^f^^^ Cap) = we can 

simply substitute 

without changing the validity of (^5|). Hence it follows that 
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+ 2ii^ - ») - "';ii-;p"' - * * + °<^"''' 

in agreement with (p^). 

To summarize, up to the Darwin order the expansion in powers A^/^ of the 
Vlasov-Maxwell system does agree with the infinite particle number hmit of the 
particle system's effective dynamics. 



3.2 Approximation to the order of radiation reaction 

The purpose of this section is to push the expansion of the Lorentz force one 
step further to include the terms of order 0{\^^^), i.e., we consider 

FL^Ea + X^/'^Ei + \{E2 + u x Bi) + \^''^{E:i + v x B2) + 0{X^), (37) 

cf. (p7|). The equations to be solved for the coefhcients E^ and B2 are divi?3 = 
47rp3, curli^s = —dtB2, divB2 — 0, and curli32 — dtEi + Anji, and this leads to 

AE3 = 47rVp3 + dt {dtEi + Anji ) , AB2 = -47r curlji . 

The solution of the latter equation is analogous to that of ( p^ and just con- 
tributes to the expansion of j. Concerning E3, it could be treated just as E2 
was, obtaining a solution which vanished at infinity. There would be no trace 
of radiation reaction and this solution is not the appropriate one for obtaining 
an approximation of a solution of the Maxwell equations at the 1.5PN level. 
The equations of the pure post-Newtonian expansion of section 2 do not con- 
tain information about damping in this case. They need to be supplemented by 
asymptotic conditions different from the condition of vanishing at infinity. In the 
case of the Maxwell equations this additional information could be obtained by 
doing expansions of the fundamental solution different from the post- Newtonian 
one. 



As in the previous section 3.1, the connection can be made between this 
result and the N 00 limit of a system of N particles. In ||l^ it was proved 
that to the order of radiation reaction the particle dynamics is governed over 
long times by 

^TkinK) = GSD(g,«,«) -f G^R(g,z;,i;) + 0{\'), 1 < a < N, 
cf. (H), with the new C'(A3/2)-term 

- ^^^^7^'^-^^^'^ (..-..0) (38) 
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accounting for the radiation reaction. We do assume now that there are two 
species of particles of equal unit mass and of total number N — 2M, such that 
e„ = e = 1 for a = 1, . . . , M and Ca = -e = -1 for a = A/ + 1, ... , 2M. The 
respective particle number densities are denoted by and Decomposing 

N 

/3,/3' = l,/3#/3' 



the X]fl^/3'=i B^B' ™ accordingly, we see that 



-^X^ J J dxidpi J J dx2dp2 f^{xi,pi)f {x2,P2) 

1 f X Hxi - X2) ■ {pi - P2) , A 

-(P1--P2) . 7^ {X1-X2)]. (39) 



\xi - X2\^ \xi - a;2p 

On the other hand, recalling the dipole moment D{t) from (H) and D^'^\t) from 
T3I) we have 

[ [ p{f+ - r)dxdp, D{t)^D^'\t) + 0{e). 



Differentiating the latter relation w.r.t. time, and using the explicit form of 
the electric field VU obtained from AU — Airp in conjunction with the Vlasov 
equations ( p] ) and ( [l^ ) for /+ and we obtain through explicit evaluation 
that 

Xl ~ X2 



1 1 1 1 c^Pi dx2 dp2 {p2 ■ Va;^! 



\xi - X2\^ 

dxi dpi (pi • V)VC/ (/+ + /-) + 0(e) 
dxi dpi dx2 dp2 ( -. pripi - P2) 

\\XI - X2\'^ 

2>{xi - X2) ■ {pi ~P2) , . 
\ \i 2^1 - X2) 

•(/+(X2,P2) - r (X2,P2)) + r + 

Now we observe that using the change of variables x\ <-> X2 and p\ <-> p2 it is 
clear that the term containing f^{x2,P2)f^ixi,pi) as well as the one containing 
f~{x2,P2)f~ixi,pi) does vanish, whereas the other two terms arc the same and 
add up. Consequently we finally find 

D {t) = ^ J J J J dxi dpi dx2 dp2 ^ 1^^ ^^^|3 (Pi ~P2) 

3(xi -X2) ■ jpi -P2) I A . . 

^ {Xx~ X2)\r(xx,px)j [X2,P2)- 

Comparing this result to ( p9[ ) we see that the VPD model considered in section 
H is a quite reasonable one to describe electromagnetic radiation reaction. 
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4 Outlook 



In this paper various models of radiation damping and their mutual relations 
have been discussed. Potential difficulties in establishing mathematical results 
about these models have been pointed out. In the models of section 1 the force 
on a particle is determined instantaneously by the distribution of all particles. 
These were self-contained models of continuum mechanics. Another type of 
model for radiation damping has been studied in the literature. In that case 
the description of matter is a schematic one (one or more oscillators) while the 
field satisfies a wave equation. Examples are a model of Burke |^ and one of 
Aichelburg and Beig whose mathematical properties were further studied in 
and It would be interesting to make connections between the latter 

models and those discussed in the present paper. 

The central open question concerning the models for electromagnetic radi- 
ation damping in section 2 is that of proving that solutions of these models 
approximate solutions of the Vlasov-Maxwell system in some appropriate sense. 
The results of section 3 may help to give indications how results of this kind 
should be formulated and proved. However all that was done up to now was 
to point out a formal similarity between the equations. If a complete rigorous 
analysis of the case of the Vlasov-Maxwell system (or even of a model such as 
that of |l) were obtained then it could serve as a basis for an attack on the 
much harder case of radiation damping for the Einstein equations. 

Acknowledgements: The authors are grateful to H. Spohn, M. Kiessling and 
B. Schmidt for discussions. 
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